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The DO collaboration reports direct evidence for electroweak production of single top quarks 
through the i-channel exchange of a virtual W boson. This is the first analysis to isolate an 
individual single top quark production channel. We select events containing an isolated electron 
or muon, missing transverse energy, and two, three or four jets from 2.3 fb _1 of pp collisions at 
the Fermilab Tevatron Collider. One or two of the jets are identified as containing a b hadron. We 
combine three multivariate techniques optimized for the i-channel process to measure the t- and 
s-channel cross sections simultaneously. We measure cross sections of 3.14+q' 94 , pb for the t-channel 
and 1.05 ± 0.81 pb for the s-channel. The measured t-channel result is found to have a significance 
of 4.8 standard deviations and is consistent with the standard model prediction. 

PACS numbers: 14.65.Ha, 12.15.Ji, 13.85.Qk 



The DO and CDF collaborations at the Fermilab 
Tevatron pp Collider have recently observed electroweak 
production of single top quarks [H,0] > measuring the total 
single top quark production cross section as well as the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix element Q 
|Vtb| directly. In the standard model (SM), the two 



main production modes at the Tevatron resulting in 
a single top (or antitop) quark final state are the t- 
channel exchange of a W boson shown in Fig. QJt and 
the s-channel production via the decay of a virtual 
W boson. The two observation analyses measured only 
the combined single top quark cross section, assuming 
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the SM ratio of the two production modes. This ratio 
is modified in several new physics scenarios, for example 
in models with additional quark generations, new heavy 
bosons [1, H, H, 0| , flavor-changing neutral curre nts I j, M , 
Hoj . or anomalous top quark couplings [Til. [l2l, [i~3l. fl4Ul5|| . 
In this Letter we remove this assumption and use the 
t-channel characteristics to measure the t-channel and 
s-channel cross sections simultaneously, thus providing 
a t-channel measurement independent of the s-channcl 
cross section model. The main characteristic of the t- 
channel which separates it both from the s-channel and 
the backgrounds is the pseudorapidity distribution of 
the light quark jet, shown in Fig. [T]d. The predicted 
cross section for t-channel (s-channel) production is 
2.34 ± 0.13 pb (LI 2 ± 0.04 pb) for a top quark mass 
m = 170 GcV [3. 
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FIG. 1: Representative Feynman diagram for t-channel 
single top quark production and decay (a) and parton-level 
pseudorapidity distribution of the final state objects in top 
production (excluding antitop), requiring each object to have 
transverse momentum > 15 GeV (b). 

This analysis extends the DO single top evidence [l7], 
liljj and observation analyses [l|, utilizing the same 
dataset, event selection, and signal/background modeling 
as the observation analysis, but training multivariate 
filters specifically to extract t-channel single top quark 
events. We use 2.3 fb _1 of data collected by the DO 
experiment [l9| at the Fermilab Tevatron pp Collider 
between 2002 and 2007 (Run II). The measurement 
selects final states containing one high transverse 
momentum (j>t) isolated lepton (electron or muon), large 
missing transverse energy (]$t), a b quark jet from the 
decay of the top quark (t— *Wb— >ivb) , a light quark 
jet produced in association with the top quark, and 
a spectator b quark jet from gluon splitting in the 
initial state. We allow for one of these jets not to be 
identified as well as for the presence of an additional jet 
from gluon radiation. The backgrounds are W bosons 
produced in association with jets, it pairs, and multijet 
production, where a jet is misreconstructed as an electron 
or a heavy-flavor quark decays to a muon that satisfies 
isolation criteria. Z+jets and diboson processes form 
minor additional background components. We treat s- 
channcl single top quark production as a background 
during the multivariate training but measure its cross 
section simultaneously with the t-channel measurement 
as explained below. 

We look for t-channel and s-channel single top quark 



production in events with two to four jets with px > 
15 GeV and pseudorapidity \rj\ < 3.4, with the leading 
jet additionally satisfying pr > 25 GeV. We require 
20 < J/jt < 200 GeV for events with two jets and 
25 < J/jt < 200 GeV for events with three or four 
jets. Events must contain only one isolated electron with 
Pt > 15 GcV and \r]\ < 1.1 (pr > 20 GcV for three- or 
four-jet events), or one isolated muon wither > 15 GeV 
and I -77 1 < 2.0. The background from jets misidentificd as 
leptons is kept to approximately 5% by requiring the total 
transverse energy of all final state objects Ht{£,$t, jets) 
to be greater than 110 to 160 GeV, depending on the 
analysis channel, and by demanding that the JpT is not 
along the direction of the lepton or the leading jet in 
the transverse plane. To enhance the signal fraction, 
one or two of the jets are required to originate from 
b hadrons, as implemented through a neural network 
(NN) 6-jet tagging algorithm [2d|. We divide the dataset 
into 24 independent analysis channels (separated by data 
taking period, lepton type, 6-tag and jet multiplicity) and 
combine the results to maximize the signal sensitivity. 
Details on the selection criteria and background modeling 
are given in Ref. [HI]. 

We generate t-channel and s-channel single top 
events with the SingleTop Monte Carlo (MC) 
generator [2l[. The kinematics of the generated events 
closely match those pre dicted by next-to-leading-order 
(NLO) calculations [22|, particularly those including 
NLO corrections to the t-channel 2 — ► 3 process shown 
in Fig. Pit t fipl . The ALPGEN leading-order MC event 
generator [24j , interfaced to pythia [25| , is used to model 
ti, VF+jets, and Z+jets background events. We use the 
CTEQ6L1 parton distribution functions (CTEQ6M for 
single top) [26[ and set the top quark mass to 170 GeV. 
We use GEANT (2?| to simulate the response of the 
DO detector to the MC events. The ti background 
is normalized to the predicted cross section [28]. The 
Z+jets contributions are normalized to NLO cross 
sections [29|. The VF+jets background normalization, 
jet flavor composition, and jet angular distributions are 
obtained from data samples. We model the background 
from multijet production where a jet is misidentificd as 
an isolated electron or muon using events from data 
containing lepton candidates which pass all of the lepton 
identification requirements except one, but otherwise 
resemble the signal events. We use pythia to model 
diboson production. 

We select 4519 lepton+jets events with at least one 
6-tagged jet, which are expected to contain 130 ± 17 t- 
channcl (93 ± 14 s-channel) single top events with an 
acceptance of (2.5 ± 0.3)% ((3.7 ± 0.5)%). The expected 
sample composition is shown in Table [T] 

Systematic uncertainties in the signal and background 
models are described in detail in Ref. [HI. The main 
uncertainties are due to the jet energy scale (JES) 
corrections and the tag-rate functions (TRF), with 
smaller contributions from MC statistics, correction for 
jet-flavor composition in LF+jets events, and from the 
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TABLE I: Number of expected and observed events in 2.3 fb~ 
for e and /i, and one and two &-tagged analysis channels 
combined, with uncertainties including both statistical and 
systematic components. The t-channel and s-channel contri- 
butions are normalized to their SM expectation. 

Source 2 jets 3 jets 4 jets 



t-channel 
s-channel 
W+jets 

Z+jets and dibosons 
tt -> U, ^±jets 
Multijets 
Total prediction 
Data 



77 ± 10 
62 ±9 
1829 ± 161 
229 ± 38 
222 ± 35 
196 ± 50 
2615 ± 192 
2579 



39 ±6 
24 ±4 
637 ±61 
85 ±17 
436 ± 66 
73 ±17 
1294 ±107 
1216 



14 ±3 
7±2 
180 ±18 
26 ±7 
484 ± 71 
30 ±6 
742 ± 80 
724 



IF±jets, multijets, and tt normalizations. The total 
uncertainty on the background is (8-16)% depending on 
the analysis channel. Uncertainties on JES, TRFs and 
the modeling of W±jets kinematics affect not only the 
normalization but also the shape of the discriminant 
distributions. Since the VF±jets background normal- 
ization and kinematics are constrained by data, theory 
prediction uncertainties for IF±jets are negligible 
(including those on the parton distribution functions and 
factorization and normalization scale). 

We apply three independent multivariate analysis 
techniques to separate the small t-channel single top 
signal from the large backgrounds, based on boosted 
decision trees (BDT) l3Cl l3ll . |32| . Bayesian neural 
networks (BNN) 133. l34| , and the matrix element 
(ME) method [35L l36|. These techniques and their 
application are described in detail in Ref. For 
this analysis we use the same set of variables as in the 
observation Letter However, only t-channel single 
top events are considered signal during the optimization, 
whereas s-channel single top events are included in the 
background, normalized to the SM expectation. Fig. [5] 
shows comparisons between the t-channel signal, the 
background model, and data for the three individual 
discriminants. 

The three multivariate techniques use the same data 
sample but are not completely correlated. Their 
combination leads to increased sensitivity and a more 
precise measurement of the cross section. We achieve 
this by training a combination BNN which uses the 
three individual discriminant outputs as inputs. Fig. [3] 
shows the combination discriminant output for data 
superimposed on the background and signal models. 

We verify the accurate modeling of the data in 
background-dominated control regions for the two main 
background categories. Fig. [4ji shows the t-channcl 
discriminant in a W+jets dominated sample of 2-jet, 1- 
tag events with Ht < 175 GeV. Fig. 0b shows the t- 
channcl discriminant in a tt dominated sample of 4-jet, 
1-tag or 2-tag events with Ht > 300 GeV. These studies 
confirm that backgrounds are well-modeled across the full 
range of the discriminant output. 

We use a Bayesian statistical analysis [37| to measure 
the production cross sections. In a first step we compute 
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FIG. 2: Comparison of the signal and background models to 
data, for the BDT discriminant (a,b), the BNN discriminant 
(c,d), and the ME discriminant (e,f), for the full discriminant 
range (a,c,e) and the signal region (b,d,f). The bins have 
been ordered by their expected t-channel signahbackground 
ratio and t-channel and s-channel single top distributions are 
normalized to the measured cross sections. 
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FIG. 3: Comparison of the signal and background models 
to data for the combination discriminant output, for the full 
range (a) and only the signal region (b). The bins have 
been ordered by their expected t-channel signahbackground 
ratio and t-channel and s-channel single top distributions are 
normalized to the measured cross sections. 
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the two-dimensional posterior probability density as a 
function of both t-channel and s-channel single top quark 
cross sections. The combination discriminants for t- 
channcl and s-channel single top, remaining background, 
and data are used to build a binned likelihood as a 
product over all analysis channels and bins. We assume 
a Poisson distribution for the observed counts, and flat 
prior probabilities for positive values of the t-channel and 
s-channel signal cross sections. Systematic uncertainties 
are described by Gaussian priors, and their correlations 
amongst all bins in all channels are preserved. The 
posterior probability density is shown in Fig. [5] Also 
shown are the SM expectation as well as several represen- 
tative new physics models to illustrate the sensitivity 
of this analysis. Dedicated searches should be able to 
address flavor-changing neutral currents with a Z boson 
coupling to the top and up quark with a strength of 4% 
of the SM coupling jij or a top-color model with a tb 
bound state (Top Pion) with a mass of = 250 GeV [1] , 
while a 4-quark-generations scenario with CKM matrix 
element |V ts | = 0.2 [38[ or a top-flavor model with new 
heavy bosons at a scale m x = 1 TeV Q will be more 
challenging to identify and might have to wait for LHC 
studies. 
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FIG. 5: Posterior probability density for t-channel and s- 
channel single top quark production in contours of equal 
probability density. Also shown are the measured cross 
section, SM expectation, and several representative new 
physics scenarios [3, [3|| ■ 



In a second step we obtain the t-channel posterior 
probability density from the two-dimensional posterior 
in Fig. [5] by integrating over the s-channel axis, thus 
not making any assumptions about the value of the s- 
channel cross section. We have analyzed ensembles of 
pseudo-datasets generated at several different t-channel 
and s-channel cross sections to verify the linearity of the 
measured t-channel cross section and its independence 



of the input s-channel cross section. From the t-channel 
posterior we extract the cross section and uncertainty 
for t-channel single top quark production as 3.14lg gQ pb. 
We similarly extract the s-channel cross section as 1.05 ± 
0.81 pb by integrating over the t-channel axis. 

We compute the significance of the t-channel cross 
section measurement using pseudo-datasets generated 
from the background model (including SM s-channcl 
single top) and taking all systematic uncertainties 
into account in a log- likelihood-ratio approach 0, [39| . 
For each pseudo-dataset we calculate the ratio of the 
probabilities for two hypotheses: that the pseudo-dataset 
is described by the background model only (including 
SM s-channel), and that it is described by SM t-channel 
single top plus backgrounds. We measure the p- value 
by counting the fraction of background-only pseudo- 
datasets with a ratio that is more signal-like than the one 
observed in data. The observed p-value is 8.0 x 10~ 7 , 
corresponding to a Gaussian significance of 4.8c, and 
the expected p-value is 9.7 x 10~ 5 , corresponding to a 
Gaussian significance of 3.7a. 
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FIG. 6: Ht (a), reconstructed top quark mass (b), light quark 
jet pseudorapidity multiplied by lepton charge (c), and t- 
channel top quark spin correlation (d, see text) for events 
with a ranked combination output > 0.91. The t-channel 
and s-channel contributions have been normalized to their 
measured cross sections. 

We have checked the consistency of the observed signal 
with SM t-channel events in several kinematic distri- 
butions. Fig. [5] shows comparisons between the observed 
data, the background model, and the t-channel signal 
for four different kinematic distributions for events with 
a ranked discriminant output > 0.91. Shown are four 
important kinematic variables for t-channel single top 
quark production: Ht] the reconstructed top quark 
mass; the lepton charge multiplied by the pseudorapidity 
of the leading non-6-tagged jet (cf. Fig. QJ>); and the t- 
channel spin correlation in the optimal basis 0, SH, i.e. 
the cosine of the angle between the light quark jet and the 
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lcpton, both in the reconstructed top quark rest frame. 
While the background shapes resemble the signal in the 
high ranked discriminant output region, the presence of 
the t-channel signal is nevertheless clearly evident in each 
distribution. 

In summary, we have presented the first direct evidence 
of the t-channcl mode of single top quark production 
using 2.3 fb" 1 of data at the DO experiment. We measure 
a t-channel cross section of 3.14+g'gg pb and a s-channel 
cross section of 1.05 ± 0.81 pb. The measured cross 
sections are consistent with the SM expected values. The 
observed i-channel signal corresponds to an excess over 
the predicted background with a significance of 4.8 a. 
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